Spin Dependence of Interfacial Reflection Phase Shift at Cu/Co Interface 
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The spin dependent reflection at the interface is the key element to understand the spin transport. 
By completely solving the scattering problem based on first principles method, we obtained the spin 
resolved reflectivity spectra. The comparison of our theoretical results with experiment is good in 
a large energy scale from Fermi level to energy above vacuum level. It is found that interfacial 
distortion is crucial for understanding the spin dependence of the phase gain at the Cu|Co interface. 
Near the Fermi level, image state plays an important role to the phase accumulation in the copper 
film. 
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I. INTRODUCTION 

The spin dependent transparency of the nonmagnetic 
metal|ferromagnet (NM|FM) interface plays a crucial role 
in the magnetoelectronics. The spin dependent trans- 
mission and reflection result in the spin dependent in- 
terface resistance, which dominates the giant magneto- 
resistance effect 1 -. More interestingly, the spin dependent 
reflection at the NM|FM interface also yields the spin 
torques 2 ' 3 at the interface. Spin torque across NM|FM 
can be formulated in terms of the mixing conductance 3 
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where M is the total number of channels (eigenstates) , 
and \exp(i<f>™n\) denotes the reflection from chan- 
nel m to channel n. The reflection phase gain difference 
(PGD) defined as <fi™ n — 4>" m plays an important role to 
understand the spin torque. 

However, measuring spin dependent transport through 
FM|NM is far from trivial, since the electrons with differ- 
ent spin orientation flow together in most of the electric 
transport process. It is possible to probe the spin de- 
pendent transmission by making use of the two current 
model^. The Andreev reflection (AR) spectroscopy at 
ferromagnet-superconductor interfaces in ballistic point 
contacts (PCAR)&£, as an alternative way to measure 
the spin polarization, contains all the information of the 
scattering matrix at the interface but experimental data 
are very difficult to be resolved^. It is deserved to have 
one method which can measure the transmission or re- 
flection through an interface for different spins separately. 

Recently, Wu et al.^ have performed the measure- 
ment of the spin-dependent electron reflection on the 
Vacuum(Vac)|Cu|Co(001) film. Spin polarized incident 
electrons were employed to detect spin dependent reflec- 
tion at Cu|Co interface. The PGD between two spin 
channels at Cu|Co interface was obtained by measuring 
the relative shift of reflectivity spectra for different spin 
polarization direction. Such an experiment provides us 
an opportunity to test the theoretical calculation of elec- 
tronic transport with experiment in great details, espe- 



cially the subtle reflection phase. Furthermore, the effect 
of atomic structure of interface on the spin transport can 
be addressed. 

In this paper, the phase accumulation and PGD in 
Vac|Cu|Co(001) system is obtained from first principles 
calculation and compared with experiment in a large en- 
ergy scale. The interfacial distortion is found to be im- 
portant for understanding the spin dependence of phase 
gain at Cu|Co interface. Effects due to the distortion on 
the spin transport are also discussed. Besides, we also 
estimate the possible phase gains from image states at 
the surface of Cu film. Such extra phase gain could be 
obvious near Fermi level. 

This paper is organized as following: In Sec. II, we give 
a brief introduction of our theoretical method. The mod- 
elling of interfacial relaxation and the optimized struc- 
ture are specified. The experiment details will also be 
given in this section. In Sec. Ill, we compare the theo- 
retical results with experimental data, where interfacial 
alloy and image potential are taken into account. We 
also discuss the effect of structure relaxation on the spin 
transport at interface. Our results will be summarized in 
Sec.IV. 



II. THEORETICAL METHOD AND 
EXPERIMENT 

Several works — ili tried to calculate the phase gain 
at the NM|FM interface, but the definition of the phase 
gain at a single interface for a given energy is no unified 
due to the arbitrary phase when solving the Schrodinger 
equation. Wei and Choi*i£ obtained the phase gain at 
Ag|Fe interface by comparing quantum well (QW) states 
positions of free standing Ag film and Ag film on top of 
Fe. While An et al.r^- defined the phase gain at Cu|Co 
interface by a special fitting wave function. In fact, what 
can be measured in experiment is not the phase gain at 
single interface, but the sum of the phase gains at two 
opposite interfaces or the phase gain difference of two 
spin channels at the same interface. 
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Based on a recently developed first principles 
method- 2 , we obtained the phase gains by calculating 
the scattering matrix at interfaces. The arbitrary phase 
can be removed by noticing that the incoming wave of 
one interface is the outgoing wave of another interface 
for electrons scattering between two opposite interfaces. 
The phase accumulation obtained by the present method 
can be used to reproduce the positions of QW states from 
the electronic structure. 

In this method, the bulk and the interface potentials 
are determined in the framework of the tight binding 
(TB) linearized muffin tin orbital (MTO) with surface 
green function (SGF) based on density functional the- 
ory (DFT) and local density approximation(LDA). Dis- 
ordered systems^ are treated using the layered coherent 
potential approximation (CPA). As SGF is incorporated 
into TB-LMTO, the electronic structure of interfaces and 
other layered system can be treated in our framework. 
The potentials^ obtained in this way were used as in- 
put to a TB-MTO wave-function-matching (WFM) cal- 
culation of the transmission and reflection coefficients be- 
tween Bloch states on either side of the interface. The 
details of WFM can be found in Ref.[I2|. 

The interfacial distortion was calculated based on 
Vienna ab. initio simulation package (VASP) 15 with 
ultrasoft pseudopotentialsi^ and generalized gradient 
approximation^ for exchange correlation energy. Com- 
mon lateral lattice constant used is the average value of 
the lattice constant of both metals obtained in VASP, 
viz. a cu = a co = 3.595A. The distortion was limited 
along [001] direction. For Cu|Co interface, the modelling 
is based on the supercell structure which consists of 6ML 
fee Co(001) and 6ML Cu(001) with 12Avacuum inserted 
between them. The three layers closed to vacuum at Co 
side and the three layers closed to vacuum at Cu side are 
fixed. The in-plane lattice constant is fixed to the com- 
mon lattice constant. Energy cut-off is set to 400 eV and 
k grid is 19x19x1. For a given supercell the structure 
is relaxed to reach the minimum energy. The cell size is 
then varied, and a curve of structure vs energy can be 
obtained. The optimized structure with the lowest to- 
tal energy is deduced from this curve. For the surface 
Vac|Cu, the supercell structure consists of 6ML Cu(001) 
and 12Avacuum. Only the three Cu layers closed to vac- 
uum are allowed to move. The optimized structure of 
Cu surface is obtained when the supercell reaches it min- 
imum energy after relaxation. 

Due to work function mismatch, the distance between 
Cu and Co layers is stretched to 1.86A. At Cu side the 
distance between the neighboring layers is stretched by 
3%. At Co side, the distance between the topmost two 
layers is contracted by 5% and the distance between the 
second and third layers is contracted by 1%. At Cu sur- 
face, the distance between the topmost two surface layers 
is contracted by 1%. With the atomic structure obtained 
from VASP as the input into the calculations of elec- 
tronic structure and transport. The volume of atomic 
sphere at interface is a little different from that in bulk 



and selected to satisfy the principle of full space filling. 
The unrelaxed part of the two leads of interface are con- 
catenated with bulk, where electronic structure of bulk 
has been obtained independently and the self-consistent 
potential of interface are treated as the " embedding po- 
tential" of the bulk system with aid of SGF—. We found 
the relaxation of Co side is crucial to understand the spin 
dependent electron reflection correctly. 

The spin dependent normal reflection with normal in- 
cidence in Vac|Cu|Co(001) system is calculated. As there 
is only one channel of Ai symmetry in copper at the T 
point in the 2-dimcnsional Brillouin zone (2D BZ), the 
total reflectivity back to the vacuum for each spin is 

_ r% + er 4d l x r 2 c - 2e- 2d ' x r B r c cos{2k e d - <j) B - <j>c) 
~ 1 + e- Ad / x r 2 B r 2 c - 2e- 2d / x r B r c cos{2k e d - (j) B - 4>c) 

(1) 

where k e = k B z — k, k B z is the BZ wave vector, k is the 
electron momentum vector. r B exp(i(j) B ) and rcexp(i(f>c) 
are the reflection coefficients from Cu side at the Cu|Vac 
and Cu|Co, and d is the thickness of Cu film. Due to 
the inelastic scattering and the impurity scattering in 
Cu film, a mean free path (MFP) A has to be incorpo- 
rated, where MFP in the energy scale is assumed to be 
uniform to minimize our adjustable parameter in spite 
of possible energy dependence. By this way, only elastic 
reflection are taken into account and no inelastic reflec- 
tion will be considered. The reflectivity asymmetry is 
defined as (R p — R ap )/(R P + R ap ) with R p {R ap ) denot- 
ing total reflectivity with incident spin parallel to major- 
ity(minority) spin in Co. 

To justify the theoretical calculation, we made the 
comparison between our calculation and the experimen- 
tal data quantitatively. The experimental electron reflec- 
tion spectrum for the electron energy higher than the vac- 
uum level was obtained by the spin-polarized low-energy 
electron microscopy (SPLEEM) measurement. The ex- 
periment was performed at Lawrence Berkeley national 
laboratory. The detailed experimental description can be 
found in Ref.[9j|. In the SPLEEM, a spin polarized elec- 
tron beam is directed at the sample surface at normal in- 
cidence and the spin-dependent electron reflectivity can 
be measured simultaneously during the Cu growth on 
5ML Co|Cu(001). As shown in Eq.Q, the maximum 
electron reflectivity takes place at the interference condi- 
tion of 2k e d — <f) B — cf>c — 2irn, where n=integer. So by 
measuring the thickness dependent electron reflectivity 
at certain electron energy, the k e and the total phase gain 
4>b + 4>c can be derived experimentally 18 . However, the 
phase gain calculated using this method has larger error 
bar, which is larger than the PGD. To obtain the PGD 
precisely, we can determine the peak position difference 
Ad in the reflection spectrum of spin-up and spin-down 
electron, then the PGD can be calculated by the formula 
Acj) = 2k e Ad. 

The QW states for minority spin below the Fermi sur- 
face were measured by the angular resolved Photoemis- 
sion at beamline 7.0.1.2 of the Advanced Light Source at 
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FIG. 1: (color online) The calculated spin dependent electron 
reflectivity spectra with (a) unrelaxed interface, (b) relaxed 
interface and (c) the experimental spectra for the system of 
Vac|Cu(7ML)|Co(001). (d) The calculated reflection asym- 
metry for the spectra shown in (a)-(c). The experimental 
spectra have been multiplied by a factor 3 for better compar- 
ison. The mean free path used in the calculation is 7ML. 



FIG. 2: (color online) (a) Total reflection phase cf>B + 4>c of 
the unrelaxed interfaces of CujCo and Cu|Vac. (b) Total re- 
flection phase 4>b + <f>c of the relaxed interfaces of Cu|Co and 
Cu|Vac. (c) Comparison of PGD of unrelaxed, relaxed inter- 
face with experimental data, (d) PGD of the two artificial 
interfaces discussed in the text. 



the Lawrence Berkeley National Laboratory, and the ex- 
perimental details can be found in Ref-pH and Ref . [l9| . 
The electron density oscillates with the Cu thickness, the 
periods will be determined by the fc e , and the phase can 
be obtained by measuring the peak position. Due to the 
high quality of our sample, the Cu thickness can be de- 
termined precisely, so the phase gain for minority spin 
can be obtained with high accuracy. 



III. NUMERICAL RESULTS AND DISCUSSION 

To identify the effect of the distorted lattice, we cal- 
culated the energy dependence of reflectivity and reflec- 
tivity asymmetry for Vac|Cu(7ML)|Co(001) system with 
unrelaxed and relaxed interfaces shown in Fig[TJ where 
the unrelaxed interface is assumed to be of idea fee lat- 
tice. The amplitude of theoretical reflectivity is not very 
sensitive to the interfacial distortion. The amplitude ob- 
tained in experiment is almost 3 factor less than theoret- 
ical results. Due to the introduction of MFP in the Cu 
film, the amplitude is dominated by the Vac|Cu surface. 
As the measurement of reflectivity is during the growth 
of sample, the discrepancy mainly results from the im- 
purities and inelastic scattering at Cu surface. However, 
the asymmetry spectra is more sensitive to the interfacial 
distortion, since only the electron reflection at Cu|Co in- 



terface is spin dependent. By quantitatively comparing 
our ab. initio results and experimental data as shown 
in FigQ] (c), we can draw a conclusion that the interfa- 
cial distortion is very crucial to obtain the correct spin- 
dependent electron reflection spectrum. Here, we found 
that A=7ML±1ML in Cu from 8cV to 16eV gives the 
best reproduction of experimental data. 

The effect due to the distortion of interface is more 
obvious for the phase gain at interface. In (a) and (b) of 
Fig. [2J the summation of the reflection phases of the two 
interfaces are shown for unrelaxed structure and relaxed 
structure respectively, where the shape of the spectra for 
two spins are similar to each other and the odd points of 
reflection phases appear at the band edge for each spin. 
The PGD spectra for the unrelaxed and relaxed Cu|Co 
interfaces are shown in Fig. [3 (c) with experimental data. 
Sign of PGD have been reversed from positive for the 
unrelaxed interface to negative for the relaxed interface, 
which is consistent with experimental results. Note that 
LDA fails to produce the image-like potential out of Cu 
surface 20 . With referring to the empirical model^i, we 
have estimated the effect due to the image states on the 
reflection phase at VacjCu by shifting the potential at 
vacuum side to reproduce the surface potential 22 . For 
energy far higher (> 8eV) than Fermi level, the effects 
due to the image states are small. 

To clarify the source of PGD for the relaxed inter- 
face, we have artificially constructed two interfaces (i) 
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TABLE I: Interface resistance SR in units of fQm 2 , the mix- 
ing conductance G m ix in units of 10 15 fi -1 m -2 . 



system 


R 1 


Rl 


Re(Gmix) 


lTYl(Grnix 


unrelaxed CuCo (001) 


0.303 


1.91" 


0.556 


-0.021 


relaxed CuCo (001) 


0.304 


2.42 


0.547 


-0.021 


relaxed CU50C050 (001) 


0.306 


1.53 


0.531 


-0.026 


unrelaxed AgFe (001) 


1.00 


6.79 


0.445 


-0.003 


relaxed AgFe (001) 


1.04 


6.10 


0.444 


-0.0006 


relaxed Ag 50 Fe 5 o (001) 


1.02 


2.59 


0.441 


0.0002 



"The slight difference of SR between present results and those in 
Ref . 1291 comes from different atomic potentials. Therein, spd basis 
is used. 



FIG. 3: (color online) The calculated reflection phase <pB+4>c 
of the majority spin(solid circle) and minority spin(open tri- 
angle) with the experimentally measured phase(open square). 
The solid line indicates the lower limit (position of error bar) 
of minority spin after taking interfacial alloy into account as 
discussed in text. 



and (ii). In (i), the structure at Co side of the relaxed 
Cu|Co interface was replaced by the ideal fee lattice. In 
(ii), the structure at Cu side of the relaxed Cu|Co in- 
terface was replaced by the idea fee lattice. As shown 
in of Fig. fj] (d), PGD of (i) is very similar to the spec- 
trum of the completely relaxed interface, which indicates 
that Co side of the relaxed structure is the decisive fac- 
tor for the PGD of the relaxed Cu|Co interface. The 
penetration of wave into the Co side will induce the spin 
dependent phase accumulation and the distortion of this 
part could affect greatly the PGD. The main difference 
of PGD between experimental data and our theoretical 
results is the general shape of the energy dependence. 
Two odd points exist in the theoretical results, however, 
there is no such clear point in experimental data at cor- 
responding energy. Such difference can not be attributed 
to the uncertainty of interfacial structure model, as the 
odd points are due to the band edges of the two spin 
species(see Fig. [2] (a)&(b)) and the interfacial structure 
couldn't smear the mismatch of band structures of Cu 
and Co. The possible contribution to the difference may 
come from the electron-phonon (ep) coupling at interface 
which is not incorporated in our calculation yet. The 
ep coupling has been found could be important for the 
phase shift at metallic interface 2 ^. With the including 
of phonon, the energy dependence of reflection phase is 
dominated not only by the electron band structure but 
also by the phonon spectrum. It is expected that the dif- 
ference of the energy dependence of PGD may be com- 
pensated by the mixing of electron band structure and 
phonon spectrum. 

In experiment, one monolayer alloy could be located at 
interface. It was found that variation of the PGD due to 
1ML substitutional alloy Cu^Coi-^O < x < 1) inserted 
between Cu and Co is small(< 2 DEG), which is con- 
sistent with above discussion as 1 ML alloy only slightly 



changes the structure at Co side. However, interfacial 
alloy might induce some extra phase gain when electrons 
are reflected from the interface, which could change both 
energy and Cu film thickness dependent spectra of reflec- 
tivity. 

Up to now, the phase discussed is above Fermi level. 
The phase near Fermi level may be more meaningful as 
most electronic transport experiments are carried out 
around this energy. Experimentally, it is very hard to 
separate the contributions of majority and minority spin 
at Fermi level. However, it is well known that only the 
electrons with minority spin can form QW states at Fermi 
level in Cu|Co(001) system. So the phase gain of minor- 
ity spin can be obtained experimentally. Fig. [3] shows 
our theoretical calculation 28 and the experimental data. 
As we discussed, the image states may have the effect 
on the phase gain for the electron lower than Fermi level. 
The including of image states shifts down the total phase 
up to 50 DEG in this energy scale as shown by the un- 
certainty bar. Such large correction is important when 
comparing with experimental data. Our theoretical re- 
sults give the good reproduction of the slop of total phase. 
The small discrepancy in magnitude can be roughly com- 
pensated by replacing the layer of Cu at interface with 
1ML CU0.3C00.7 alloy, which could further shift down the 
phase spectra as shown by the solid line in Fig. [3l The 
total phase of majority spin is also present, which is far 
from the experimentally measured one. 

For r point in 2D BZ, we have shown that the relax- 
ation of atomic structure could greatly affect the spin 
transport. But, this is not enough to understand the 
spin transport. To discuss the effect of interfacial dis- 
tortion on the quantities such as interface resistance or 
mixing conductance, we should extend our calculation to 
the whole 2D BZ at Fermi level. The results are shown 
in Tab. HI For Cu|Co interface, interface resistance of 
majority spin is less affected, however, minority spin is 
well affected by structure relaxation. For majority spin, 
d bands in Cu and Co are fully filled and those states at 
Fermi level are from s bands which are essentially free 
electron like states. The interface could be equivalent to 
a very low potential step for those states which have the 
same symmetry in both materials. Interfacial distortion 
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is just like the variation of the potential profile, which is 
expected to have little effect on the transport. For minor- 
ity spin, s band still dominates in Cu, but d bands in Co 
are only partly filled, which lead to the multiple sheets 
of Fermi surface in Co. Interface resistance should be 
attributed to the mismatch of band structure in Cu and 
Co^, e.g. mismatch of Fermi velocity and mismatch of 
symmetry of wavefunctior>22. Interfacial distortion could 
play important roles in those issues. Furthermore, our 
results show that one layer of interfacial alloy could en- 
hance such effect. This can be understood as more states 
of minority spin at Co side could be diffusively coupled 
to those at Cu side with aid of the alloyi^. 

However, G m i X is not sensitive to interfacial distortion. 
According to the definition of G m i X , there exists a term 
r" m \\r" m \exp(i(4>^ rn — <^" m )) with the summation of the 
2D BZ. For the single point in 2D BZ this term is sensi- 
tive to the interfacial distortion, however, the summation 
could lead to the cancellation of this term among differ- 
ent points in 2D BZ due to the phase factor (f>™ m — 0™ m . 
The real part of G m i X obtained in our calculation is very 

2 

close to the value of %-M, which means strong cancel- 
lation did happen. As a result, the sensitive issue of 
interface reflection phase could be trivial one. 

To confirm above argument, we also studied the 
Ag|Fe(001) interface. The modelling of the interfacial 
relaxation is just like that for Cu|Co interface, but note 
that in this case we need to rotate the bee Fe(OOl) plane 
by 45 DEC to match the fee Ag(OOl) plane. The lattice 

constant ap e — ^aA g — 2.855A. After relaxation the 
distance between Fe and Ag layers is about 1.39 times 
that of neighboring layers' distance in bulk Fe. The con- 



traction of the top two layers at Fe side is by 1% and 
that at Ag side is less than 1%. Similar results were 
also obtained in this system. The interface resistance of 
majority spin and the real part of G ralI are both less sen- 
sitive to the interfacial distortion shown in Tab. |TJ The 
magnitude of the imaginary part of G m ix in this case 
greatly decreases for the relaxed interface, but the imag- 
inary part is too small compared with the real part to 
contribute to spin transfer torque 3 . 



IV. CONCLUSION 

The phase information obtained from experiments can 
be well understood by our ab. initio calculation in a large 
energy scale, it was found that structural distortion at 
interface is important for theoretical model to interpret 
experiment. The PCD of two spins obtained in our calcu- 
lation is comparable with experimental results. At Fermi 
level, extra phase gain due to image states becomes ob- 
vious. After integration of 2D BZ, we also found that 
the interface resistance of minority spin can be well af- 
fected by the interfacial distortion, however, the interface 
resistance of majority spin and the PGD related quan- 
tity mixing conductance are not sensitive to the interface 
structure. 
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